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There has been a recent tendency to apply Schrodinger's wave equation to macro- 
scopic domains, from Bose-Einstein condensates in neutron stars to planetary or- 
bits. In these applications a hydrodynamical interpretation, involving vortices in 
some 'fluid' medium is often given. The vortex picture appears surprising in light of 
more traditional interpretations of Schrodinger's equation, and indeed often appears 
to rely on ad hoc analogies. The purpose of this letter is to examine the vortex 
hypothesis in light of a simple, transparent mathematical framework. We find that 
Schrodinger's equation implies waves can be vortices also for a class of wave functions. 
Vortices occur in pairs that perform a sort of quantum computation by collapsing 
into either or 1. Vortices collapsing to ('0- vortices') are longer-lived, and their 
ratio to 1- vortices is discussed. 
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Introduction 

Can a wave be a vortex? The circular motion of vortices is seen in spiral galaxies, 
hurricanes, and quantum fluids. Are all these, in some sense, the same thing? Is it legitmate 
to apply Schrodinger's wave equation to all of them? Is it acceptable to establish an analogy 
between Schrodinger's equation and a hydrodynamic or similar equation in order to invoke 
vortices, and then to return to Schrodinger's equation to invoke quantized observables of 
some type — whether energy levels or planetary orbits? The Bohr-Sommerfeld model of 
quantized energy levels of electrons around the nucleus, or quantized angular momentum, 
could be visualized as little planets revolving around a nucleus-sun. Is it now legitimate to 
invoke quantum theory to describe planetary orbits? Wasn't Bohr's theory superseded by 
the new quantum theory of Heisenberg and Schrodinger? 

This letter addresses a single primitive aspect of these various questions: Can 
Schrodinger's wave equation imply vortices without ad hoc appeals to other areas of physics? 
If so, the goal is to demonstrate this in a simple, transparent methematical framework. We 
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want to know if a vortex can arise from a wave in Schrodinger's wave equation. 



The wave function as the complex power z c 

Consider the wave function ip = z c , where z > is real, and c = x + iy is complex. We 
may expand ip, a complex number ip = u + iv, as follows: 

^ = u + iv = z c = e clnz = e xln z+iy ln z = z x [cos{y \nz)+i sin{y In z)\. (1) 

For ip to be analytic, we need it to satisfy the Cauchy-Riemann equations: 

(2) 



du dv 

dx dy 

du dv 

dy dx 



(3) 



where 



u = z x cos(y ln z) (4) 
v = z x sin(y ln z). (5) 



We have 



du 

— = z x ln z cos(y ln z) (6) 
dx 

dv 

— — z x ln z cos(y ln z) (7) 
du 

— = —z x ln z sin(y ln z) (8) 
dv 

— = z x ln z sin(y ln z). (9) 
dx 

Hence the Cauchy-Riemann conditions are satisfied, and we may take the complex derivative 

dtp du .dv . r . , 

That is, the derivative of ip with respect to c is equal to ip multiplied by ln z. Taking the 
derivative simply contracts or stretches the length of ip according to whether the absolute 
value of ln z is less than or greater than 1 (and changes the direction of ip if ln z is negative). 
Simularly we have 



~d& 



= (In z) 1 i> (11) 
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The second derivative of ip with respect to c simply multiplies ip by (In z) 2 . That is, (In z) 2 
is the eigenvalue of the differential operator with respect to the eigenvector ip. 
Because ip — u + iv is analytic, u and v satisfy Laplace's equation 

d 2 u ^ d 2 u _ g _ d 2 v ^ d 2 v 
dx 2 dy 2 dx 2 dy 2 

Agaih, because ip(c) is analytic, we have from Cauchy's theorem that for a simple closed 
curve C that is sufficiently smooth 



ip{c)dc=0 (13) 

c 

and the additional result that at a point c = a inside C 

Ma) = — i ^-dc. (14) 
yy ' 2m Jc c- a y ' 

Note that since ip*ip = z 2x , some restrictions are necessary in order for ip to be normal- 
izable. The integral ^e 21 ln z dx converges provided z < l,x > or else z > \,x < 0. 
Otherwise we must restrict x to x E D, where D is some finite domain. 

We now let z = z(r x ,r y ,t) where r x and r y are position coordinates, and t is time. (Note 
that we have already used x and y as components of the complex c — they are not position 
coordinates.) Thus we have 

sr = £ # (16) 
^ = £ # < 17 > 



<9 2 ?/> cd 2 z c 2 , . dz . 9 c , , dz 

M = + - (18) 

f| = + - (is) 

or^ z ar^ z l or y z l or y 



Finally, let's substitute Eqs. (15-19) into Schrodinger's wave equation 

dip k 2 ,d 2 ip d 2 i 
dt 2m dr 2 dr. 



*% = -*=& + ^\ + V<r.,r,)*. (20) 

y 



We obtain 



,.dz h 2 r d 2 z d 2 z c—l,,dz. 9 . dz x9N , z TT . . . . 

*»* + SS [ ft? + w + — <fe» + ( ft-> " - r »» = a (21) 
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Note that the multiplier on potential energy U can be rewritten 

z zc* 



(22) 



c (x 2 + y 2 ) ' 

where c* = x — iy. 

Eq.(21) can be divided into real (R) and imaginary (I) parts: 

IM + % + ^t' 2 + <^ )2)| - w^) u ^ - ° <*> (23 » 

^^^^J^ir^-m (24) 

We now make some simplifications. Let x — 1 and y = 2. We have (x 2 + y 2 ) = 5. This 
yields as our real and imaginary equations 

d 2 z d 2 z 2m,z . . . , . « 

+ l&--Z3-U(r x ,r v ) = 0(R) (25) 



<9r 2 9r 2 5^ : 

m& l r .9z, 9 . <9z .9,. 2mz TW « . 1 . . 

+( ft; ))1+ M rI/(r - ,r ' )=0(/) (6) 

Note that Eqs. (25-26) are dictated by Schrodinger's wave equation. 

For a fixed potential U(r x ,r y ) = Uf, simple solutions to the real (R) equation, Eq.(25), 

are 

z = e^^^. (27) 



The vortex 



The imaginary equation, Eq.(26) is more interesting. Again we consider a fixed potential 
Uf. The term (J^) 2 + (J^) 2 is the equation for a circle whose components are and J^-. 
The radius of the circle is determined by ^[^ff + ^-Uf]. Whether the circle expands or 
contracts — thereby creating a spiral — depends on ||. A spiral is simply a vortex restricted 
to two dimensions. So, without loss of generality — since we can easily make z a function of 
three location co-ordinates — we shall refer to a spiral as a vortex. Thus we have that the 
imaginary part of the Schrodinger equation yields for the wave function ip = z c a vortex 
solution. Waves can be vortices also. Recall in this connection Eqs. (10-11), which show 
that changes in c do not materially change the nature of the underlying wave. So our use of 
specific values for x and y in c = x + iy does not change the nature of this conclusion. In 
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the spirit of 'there be dragons here' written on medieval maps, we can say 'there be vortexes 
here'. 

One simple solution to the imaginary (I) equation, Eq.(26), is 

z = e k(r x +r y )-3k^t ( 2 g) 

where 



k = ^U, (29) 
We can see the motion of the vortex by solving for u and v in 

^ = U + iv = Z° = e Hr,+ry)-3k^t+i2k(r x +r y )-6ik^t_ ^ 

Hence 

u = e k ^ +r ^- 3k2 ^ t cos(2k(r x + r y ) - Qk 2 -t) (31) 

v = e k ^ r ' +r ^- 3k2 ^ t sin(2k(r x + r y ) - 6k 2 -t). (32) 

m 



In the (u,v) plane for if;, the initial radius \J{u 2 + v 2 ) of the vortex at t = is e k ^ rx+Tv \ 

What is the meaning of (§-J 2 + (||) 2 ? We find - Z [{£-J 2 + (||) 2 ] = 2k 2 z. So the 
derivatives of z with respect to the position co-ordinates describe a circle with radius kz\/2. 
Now, since z is changing with time t, that means the radius of the circle is changing. 
Assuming r x + r y > 0, the radius shrinks to ky/2 at time t = ^ Vx ^ y where z = 1 and the 
wave function ip = z c collapses into ip = u = 1. An alternative solution to the imaginary 
equation, Eq.(26), is z = e~ k ^ r ' x+r ' y ^ 3k2 '^ t , in which case z and the radius of the circle shrink 
to a point as t — > oo and z — > 0, and the wave function ip = z c collapses into the origin. 
Thus, for an ensemble of vortex waves, their collapse produces a binary sequence of 0s and 
Is. 

We have assumed r x + r y > 0. If r x + r y < 0, we simply multiply r x + r y by -1 and apply 
the other solution. So, without loss of generality, we will assume r x + r y > 0. 

If a vortex collapses into 0, we will call it a 0-vortex; if a vortex collapses into 1, we will 
call it a 1 -vortex. The 0- vortices and 1- vortices occur as solution pairs to the imaginary (I) 
Schrodinger equation. If they are produced in equal quantities, but have different times- 
to-live, there will be an asymmetry in the observed types of vortices, and hence in their 
collapsed remainders — more 0s than Is, or vice- versa. Since 0- vortices live longer than 1- 
vortices, there are initially more 1-vortices that have collapsed into Is, and thus more Is than 
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Os. However, since 0-vortices are longer lived, we presume their number increases relative 
to 1-vortices, until eventually the number of 0-vortices collapsing to at any time is equal 
to the number of 1-vortices collapsing to 1. Thus the number of Os and Is would occur in 
equal proportions, but there would be more 0-vortices than 1-vortices. 

The wave function ip can be normalized to unity. We will calculate normalizing constants 
A Q and Ai for 0-vortices and 1-vortices, respectively. For 0-vortices, the value of the integral 

/"oo „ „ p ~2k(r x +r y ) 

e m dt = — — 
o bk z 



/•OO p *-"vy 'x i ' y ) 

r ^ dt = e -2Hr x+ r y ) 1^ e - 6k ^t dt= e __ (33) 



So a normalizing constant is A = e k( - Tx+rv ^k^6^. For 1-vortices, the value of the integral 

f ™LZ+pL /-f 1 ^ n e 2Hr x +r y ) f 2k(r x +r y ) 2k(r x +r y ) _ 1 

m m 

(34) 

So a normalizing constant is Ai = k^j6^(e 2k ( Vx+ry ) — l) - ^. From the ratio of the squared 
normalzing constants, [Ao/Ai] 2 , we expect the ratio of 0-vortices to 1-vortices to be 
]^ e Ak(r x +r y ) _ e 2k(r x +r y )^ b ecause this would equalize probabilities in the production of Os 
and Is. 

Now consider again the gradient compoents of z. For 1-vortices, both differential opera- 
tors and have the same eigenvalue k. Thus the three-dimensional map (J^-, ■§^ ) z) 
begins at the lower boundary point (k, k, 1) and forms a straight line away from the origin 
given by (kz,kz,z) for z > 1. For 0-vortices, the eigenvalues of the differential operators 
are both —k, so the same map begins at the upper boundary point (on z) (—k, —k, 1) and 
forms a straight line into the origin given by (— kz, — kz, z) for z < 1. There is a one-to-one 
mapping between these two lines given by (kz,kz,z) — > (—k/z,—kz,l/z) for z > 1. As 
z — > oo, \jz — > 0. We can think of constructing these two line seqments as follows. We 
start with a 45-degree line, a ray, in the First Quadrant of the plane ( Jf- , J^), and snip it 
at (k, k). We reflect the shorter snipped piece around the origin to (— k, —k) in the Third 
Quadrant. The two line segments form the basis for the vortex pairs. However, if we square 
the co-ordinates, and map (k 2 z 2 , k 2 z 2 , z 2 ), the two line segments are reunited into a con- 
tinuous ray in the First Quadrant begining at the origin. Thus the two vortex solutions 
are hidden in + & by a quadratic uniting of the two line segments which are in fact 
disconneted in the ( Jr-, ) plane. 
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Energy levels in the vortex 

From Eqs. (31-32) the multiplier on time t is 

6k 2 - = % (35) 
m n 

and thus the energy E is 

E = ™U,. (36) 

The total energy E is y the potential energy Uf. This value is a function of our choices for 
x and y. 

For quantized energy levels, we can replace Uf by the potential U(E). Let {Ej} be a set 
of real numbers (eigenvalues) with Eq < E\ < ... < Ej < E, where we define j by the unit 
step function A(-), namely, j + 1 = J2iZo ~ Then we define U(E) as 

U(E) = ^E ,j = (37) 

g oo j-1 

^ = T«E - E ^ -Y, E H> 0- (38) 

iZ i=0 i=0 

Note that J2£Lo X(E — Ei)Ei is simply the sum of eigenvalues less than E. Eq.(36) is replaced 
by E = Ej. Then a solution for the 1- vortices is, from Eq.(28), 

where 

k = iw u{E) - (40) 

The potential U(E) is the largest energy eigenvalue Ej permitted by E (multiplied by a 
scale factor ^), so the solution for if) = z c is quantized with angular frequency uj — 
Upon a transition from Ej^\ to Ej, the change in k is 



Ah = ^[y/Wj - y/E~J, (41) 



so that k follows the jump process 



dk = J^I^dE- (42) 



Similar reasoning applies to 0- vortices. For a fixed k, the time evolution is a smooth function 
of t, but jumps as energy AE is released or absorbed according to the step function U(E). 
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Conclusion 

This letter has demonstrated that the ordinary linear Schrodinger wave equation generates 
vortex solutions for the class of wave functions represented by ip = z c , where z > is real, and 
c = x + iy is complex. Recently the author generalized Heisenberg's uncertainty principle to 
apply to stable amplitudes s a% p js]. Heisenberg considered only the Gaussian case where a = 
2. However, for stable amplitudes with a — > 0, quantum uncertainty occurs at a macroscopic 
level. This in itself suggests that applying Schrodinger's equation(s) to macroscopic domains 
may be legitimate. Moreover, when one finds that vortices are an intrinsic part of the 
Schrodinger equation, at least in certain circumstances, then the occurrence of vortices at all 
levels of the Universe is suggestive of a fractal structure, where correspondences are obtained 
by scaling and dimensional adjustment. (The Hausdorf dimension of a stable process s a ^ 
is a.) In this case we would not need Bohr's correspondence principle of a transition from 
the quantum to the classical level: the fractal structure itself would be the correspondence 
principle. All of this, of course, must be subjected to observation and experiment. 

In this vane, Nottale's quantization of planetary motion isn't so far-fetched. Of course 
Nottale 0] has argued for a scale-relative universe, but depletes much of the force of his 
thesis by also arguing that a must equal 2. Another Gaussian interpretation, coupling the 
Schrodinger equation to the Poisson equation and leading to a lognormal galaxy distribution, 
is described in jl . 

Neto et. al. 3j] use a Schrodinger-type diffusion equation to calculate planetary mean 
distances, the results of which agree with all known planetary orbits and the asteroid belt 
in the solar system, with three empty slots. In deriving their result they substitute for % 
the much larger parameter g* = 3.6 x 10 42 J ■ s. This suggests the domain of applicability of 
Schrodinger's equation is a matter of scaling ('as below, so above'). 

Volovik 

[6] examines vortices in superconductors, superfluids, and low-density atomic 
Bose-Einstein condensates, without making any reference to the Schrodinger wave equa- 
tion. However, a nonlinear Schrodinger equation has been explicitly offered as a model of 
superfluidity This is further indication of the universal applicability of the Schrodinger 
equation. 
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